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ABSTRACT 

We use direct iV-body simulations to investigate the evolution of star clusters with 
large size-scales with the particular goal of understanding the so-called extended clus- 
ters observed in various Local Group galaxies, including M31 and NGC6822. The 
iV-body models incorporate a stellar mass function, stellar evolution and the tidal 
held of a host galaxy. We find that extended clusters can arise naturally within a weak 
tidal field provided that the tidal radius is filled at the start of the evolution. Differ- 
ences in the initial tidal filling-factor can produce marked differences in the subsequent 
evolution of clusters and the size-scales that would be observed. These differences are 
more marked than any produced by internal evolution processes linked to the prop- 
erties of cluster binary stars or the action of an intermediate-mass black hole, based 
on models performed in this work and previous work to date. Models evolved in a 
stronger tidal field show that extended clusters cannot form and evolve within the in- 
ner regions of a galaxy such as M31. Instead our results support the suggestion many 
extended clusters found in large galaxies were accreted as members of dwarf galaxies 
that were subsequently disrupted. Our results also enhance the recent suggestion that 
star clusters evolve to a common sequence in terms of their size and mass. 

Key words: stellar dynamics — methods: N-body simulations — stars: evolution — 
globular clusters: general — galaxies: star clusters 



1 INTRODUCTION 

Recent observations have unearthed a new class of extra- 
galactic objects which have been termed extended star clus- 
ters. These have shown up in the outer regions of a range 
of galaxy types: spirals M31 (Huxor et al. 2005; Mackey 
et al. 2006; Huxor et al. 2008) and M33 (Stonkute et al. 
2008; Huxor et al. 2009), the dwarf irregular NGC6822 
(Hwang et al. 2005; Huxor et al. 2010b) and the Sculp- 
tor Group dwarf elliptical Scl-dEl (Sc22: Da Costa et al. 
2009). The few extended clusters (ECs) for which high qual- 
ity colour-magnitude diagrams (CMDs) are available appear 
old and metal-poor with no evidence for multiple popula- 
tions (Mackey et al. 2006). Thus they have similar stellar 
populations to those of typical Milky Way and M31 glob- 
ular clusters (GCs). A faint extended cluster in the giant 
elliptical galaxy NGC5128 has also recently been discov- 
ered (Mouhcine et al. 2010) with properties similar to the 
Palomar GCs found in the outer regions of the Milky Way 
(van den Bergh & Mackey 2004). 

The size scales of the ECs, generally measured as the 
projected half-light radius, rh,i (where we will use the three- 
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dimensional for the related half- mass radius), are in the 
range of 10 — 35 pc which places them between regular glob- 
ular clusters and low-luminosity dwarf spheroidal galaxies 
(Huxor et al. 2005; Da Costa et al. 2009). As such these ob- 
jects raise questions as to whether they are star clusters that 
have somehow evolved to much larger radii than is typical for 
Milky Way globular clusters or if they are more closely as- 
sociated with dwarf spheroidals, perhaps even to the extent 
of possessing some kind of dark matter component. Inter- 
estingly there are extended clusters in M31 with the same 
luminosity, and at the same distance from the centre of the 
galaxy, _R gc , as regular globular clusters but with a half-light 
radius a factor of ten or more greater (Huxor et al. 2010a). 
If the extended clusters are indeed a sub-class of globular 
cluster then we must ask how clusters that appear very sim- 
ilar in many respects can obtain quite distinct structural 
properties. 

In the case of two clusters of the same mass M at the 
same -R gc but with markedly different half-mass radii we 
must consider the prior history of both clusters, pos- 
sibly traced back to having distinct initial rh values. In 
turn this may indicate different modes of cluster forma- 
tion (Elmegreen 2008; Da Costa et al. 2009). Baumgardt 
& Kroupa (2007) have shown that variations in star for- 
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mation efficiency for clusters as they form from molecular 
clouds can lead to a range of initial half-mass radii. One of 
the goals of this paper is to investigate the ramifications of 
this for the long-term structural evolution of realistic star 
cluster models. In particular we would like to know if large 
differences in the structure of clusters at a Hubble time is 
more likely a result of different structures at birth, different 
internal evolution histories (e.g. formation of binaries) or 
possibly having experienced different external environments 
in the past. Of course a combination of these effects is also 
possible. We are also interested generally in producing N- 
body models of potential extended star clusters. 

In general terms the evolution of a star cluster is dic- 
tated by internal processes, primarily two-body relaxation, 
as well as the influence of the gravitational potential of the 
host galaxy. Historically the gross properties, such as the 
dissolution time, have been described in terms of R sc , M 
and r h (Fall & Rees 1977; Gnedin & Ostriker 1997). The 
first two can be used to define a tidal radius, r t , while the 
timescale for two-body relaxation can be derived from the 
latter two quantities. Typically the timescale on which stars 
are pushed across the tidal boundary and the cluster sub- 
sequently evaporates owing to two-body relaxation can be 
estimated from M and rh as well (see Gnedin, Lee & Ostriker 
1999 for a discussion) although as shown by Gieles & Baum- 
gardt (2008) the tidal field strength is the most important 
factor in determining mass-loss rates. For clusters orbiting 
in the inner regions of a galaxy the lifetime can also be af- 
fected by gravitational shocks from the disk and/or bulge of 
the galaxy (Gnedin & Ostriker 1997). If we assume that the 
half-mass radius of a cluster scales with the tidal radius, as 
is often done (and noting that r t decreases with time), then 
two clusters of the same mass at the same R gc will have the 
same rh. In this picture any differences in rh (assuming the 
clusters are the same age) requires that the masses in fact be 
different or that the clusters have followed different evolu- 
tionary paths. It may be possible that both clusters did not 
always reside at R sc , having different galactic orbits, with 
the extreme scenario being that one cluster was accreted 
during a galaxy merger (Mackey & Gilmore 2004; Da Costa 
et al. 2009). Another possibility is that internal factors such 
as the fraction of binaries (Gao et al. 1991) or the formation 
of an intermediate-mass black hole (Kormendy & Richstone 
1995) played a role. Prior iV-body models that have docu- 
mented the effect of such factors on the evolution of a star 
cluster include Giersz & Heggie (1997), de la Fuente Marcos 
(1997), Portegies Zwart et al. (2001), Baumgardt, Makino 
& Hut (2005) and Mackey et al. (2008). 

We will be looking at the evolution of iV-body mod- 
els that start with different sizes and models that evolve 
in tidal fields of different strength. For the latter we have 
in mind M31 and NGC 6822 as examples of galaxies im- 
posing a strong and weak galactic tidal field, respectively. 
We will be interested in the maximum half-mass (and half- 
light) radii that can be obtained at the same i? gc in these 
two galaxy types. In particular, we choose i? gc = 10 kpc. 
This is motivated by the location of the ECs in NGC 6822 
which extend out to 10.8 kpc in projection from the galaxy 
centre (Hwang et al. 2005; Huxor et al 2010b) as well as the 
observation of ECs and GCs of similar luminosity/mass at 
about this position in M31. This radius could be considered 
an extreme choice for M31 which has both ECs and GCs out 



to at least 100 kpc in projection (Martin et al. 2006, Huxor 
et al. 2008, Huxor et al. 2010c), with only the innermost 
presently-known ECs residing at ~ 10 kpc; however, select- 
ing this radius for our M31 models will give an indication 
of the maximum expected influence of a strong tidal field 
on EC evolution. Note also that our models should encom- 
pass the situation for additional ECs such as the two known 
in M33, which reside in a galaxy intermediate in mass be- 
tween M31 and NGC 6822 but again at comparable radii: 
12.4 and 28.6 kpc in projection, respectively (Huxor et al. 
2009). In this sense we expect our results to be quite widely 
applicable. 

The remainder of this paper is structured as follows. 
In Section [2] we describe the setup of the iV-body models, 
including the iV-body software and how the tidal field is 
modelled. This is followed by the results in Section and 
a discussion of the implications for observations of ECs in 
Section 2] 



2 MODELS 

The model clusters in this work are evolved using the NB0DY6 
direct iV-body code (Aarseth 2003). Gravitational forces are 
computed using a fourth-order Hermite integration scheme, 
without softening, on either standard CPU or graphics pro- 
cessing unit (GPU) hardware. NB0DY6 is the sister code to 
NB0DY4 which instead interfaced with GRAPE-6 hardware 
(Makino 2002) for rapid force calculation. The codes include 
algorithms for stellar and binary evolution as described in 
Hurley (2008a, 2008b) and deal directly with the effects of 
close dynamical encounters: perturbations to binary orbits, 
collisions and mergers, formation of three- and four-body 
subsystems, exchange interactions, tidal capture and binary 
disruption. 

We focus on a set of models starting with N = 100 000 
particles (stars or binaries). Each model includes a mass 
function, stellar evolution and the tidal field of a parent 
galaxy. Single star masses are taken from the initial mass 
function (IMF) of Kroupa, Tout & Gilmore (1993) between 
the limits of 0.1 — 30 Mq. When including binary stars the 
masses of two randomly chosen single stars are combined and 
then reassigned using a mass-ratio from a uniform distribu- 
tion in order to give the primary and secondary masses. Ini- 
tial periods are drawn from the log-normal distribution de- 
termined by Duquennoy & Mayor (1991) from observations 
of local solar-type stars and the eccentricities are assumed 
to follow a thermal distribution (Heggie 1975). Metallicity 
is set at Z — 0.001 for the stars. The initial positions and 
velocities are assigned according to a either a Plummer den- 
sity profile (Plummer 1911; Aarseth, Henon & Wielen 1974) 
or a King model (King 1966: Wo = 7) with the stars in virial 
equilibrium. 

The influence of the tidal field of the parent galaxy is 
modelled by assuming that the model cluster follows a cir- 
cular orbit at distance R sc from a point-mass galaxy of mass 
M g . At this stage we are not interested in the complicating 
factors presented by more realistic three-dimensional galaxy 
potentials. The implications of this decision will be discussed 
in Section^ The relevant equations of motion are developed 
by linearizing the galactic potential in a rotating reference 
frame centered on the cluster centre-of-mass with the a;-axis 
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directed away from the galactic centre and the y-axis in the 
direction of motion (see Giersz & Heggie 1997; Vesperini & 
Heggie 1997). A tidal radius (also called the Jacobi radius: 
Gieles & Baumgardt 2008), 



ft 



= (M/3M g ) 1/3 R e 



(1) 



can then be defined corresponding to the saddle point on the 
x-axis of the effective cluster potential. Within an ./V-body 
simulation the user has the freedom to set a length-scale 
R sc with one possible choice being that the outermost stars 
of the initial density profile are scaled to sit at r t , in which 
case the cluster is said to be Roche-lobe filling (e.g. Tanikawa 
& Fukushige 2005). The Plummer profile formally extends 
to infinite radius so in practice a cut-off at a radius of ~ 
10 rh is applied to avoid rare cases of large distance. For our 
Plummer models described below this leads to r max ~ 8rh, 
where r max is the position of the outermost star, and we 
define r max /r t as the tidal-radius filling factor. The related 
ratio rtx/rt is an important quantity to describe the structure 
of star cluster models. 

In this work we model parent galaxies with two distinct 
masses: M g = 9 x 10 9 Mq to represent a dwarf galaxy such 
as NGC 6822 and M g = 9 x 10 10 Mq to model a more sub- 
stantial galaxy such as M31. As discussed above all clusters 
are set to orbit at R gc = 10 kpc in these model galaxies. 

Six distinct models are performed. These are listed in 
Table Q] Models Nl, N2 and N3 are all evolved in the 
NGC 6822-like tidal field. They are identical in all respects 
except for using R sc = 7, 14 and 21 so that they have initial 
tidal-radii filling factors of 0.33, 0.66 and 1.00, respectively 
(the corresponding r^/rt ratios are given in Table [TJ. Model 
N2b is the same as N2 except that it starts with 95 000 sin- 
gle stars and 5 000 binaries rather than 100 000 single stars. 
We then have Model Ml which is evolved in the stronger 
M31-like tidal field for comparison. All of these models start 
with a Plummer density profile whereas Model M2 starts 
with a King profile and is also evolved in the M31-like tidal 
field. Both Ml and M2 start with tidal-radii filling factors 
of 1.00 (r max matches the tidal radius). The initial mass of 
each model is M ~ 58 000 Mq which gives an initial r t of 
129 pc for models Nl, N2, N2b and N3 compared to 60 pc 
for models Ml and M2. Each model is evolved to an age 
of 20 Gyr or until 5% of the stars remain, whichever occurs 
first - the latter only happens for M2 at an age of 17.4 Gyr. 
The simulations are performed using Tesla S1070 GPUs at 
Swinburne University. 

We will also lean on the results of some previous N- 
body simulations when evaluating our results. These include 
models KlOO-OOa and KlOO-OOb of Hurley (2007) which 
both featured 100 000 single stars evolved within a stan- 
dard Galactic tidal field (Giersz & Heggie 1997): an orbital 
speed of 220 km s -1 at R gc = 8.5 kpc (with corresponding 
M g = 9 x 10 10 Mq). The two models were setup in the same 
way and the difference of note was the formation of a long- 
lived binary composed of two stellar-mass black holes (BHs) 
in KlOO-OOb which altered the central structure of the clus- 
ter compared to KlOO-OOa. Along the same lines we will use 
the results of Mackey et al. (2008) who looked at the effect 
of a population of BH-BH binaries on cluster evolution, as 
well as models including intermediate-mass black-holes (e.g. 
Gill et al. 2008). Also mentioned will be models of 30 000 
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Figure 1. Evolution of the half- mass radius, rh, and the core- 
radius, r Cj for the NGC 6822 models Nl (dotted lines), N2 (dashed 
lines) and N3 (solid lines) noting that r c is always less than r^. 
The tidal radius of model N3 is also shown. 



single stars from Hurley et al. (2004) , mainly for illustrative 
purposes. 

It is important to emphasize that unless otherwise spec- 
ified the radii quoted will be based on three-dimensional 
data. The core radius, r c , comes from a density- weighted cal- 
culation (Casertano & Hut 1985) that is traditionally used 
in iV-body models and is not comparable to the quantity 
derived by observational methods. This has been discussed 
in the past (e.g. Wilkinson et al. 2003; Hurley 2007). For 
our purposes this is fine as we use r c as an indicator of 
the cluster dynamical state, in particular to determine if 
core-collapse has been reached, rather than to compare to 
observed results for actual clusters. We also use as the 
three-dimensional half-mass radius. Here we are mainly in- 
terested in the relative values between models. However, we 
will also provide the half-light radius, rh,i, calculated from a 
two-dimensional projection of the data, to give a reference 
point for comparing the size of the model clusters to real 
clusters. 



3 RESULTS 

The different initial filling factors of models Nl, N2 and N3 
lead to different initial half-mass radii (see Table [TJ and 
allow us to investigate the effect this has on the long-term 
evolution of star clusters residing within a weak tidal field. 
It is well established that the evolution of a star cluster 
is intricately linked to the two-body relaxation timescale 
which is typically characterized by the value at the half- 
mass radius, 



Uh oc 



N 



log0.4iV \ M 



1/2 



(2) 



4 J. R. Hurley & A. D. Mackey 



Table 1. Key information for the initial models of the new simulations performed in this work and 
also the models KlOO-OOa and KlOO-OOb from Hurley (2007). All models started with N = 100 000 
stars and M ~ 58 000 Mq . Column 1 gives the label assigned to each model where N represents the 
NGC6822-like tidal field and M the M31-like tidal field. In Columns 2 and 3 we show the density 
profile and the primordial binary fraction, respectively. Columns 4-8 give the half-mass radius, 
tidal radius, half-mass to tidal radius ratio, tidal radius filling factor and half-mass relaxation 
timescale at T = OGyr. All lengths are in pc units and relaxation times are in Myr. 



Label 


profile 


fb 


<h 


rt 




fmax / Tt 




Nl 


Plum 


0.00 


5.4 


129.0 


0.042 


0.33 


980 


N2 


Plum 


0.00 


10.9 


129.0 


0.084 


0.66 


2770 


N2b 


Plum 


0.05 


10.9 


129.0 


0.084 


0.66 


2697 


N3 


Plum 


0.00 


16.3 


129.0 


0.126 


1.00 


5090 


Ml 


Plum 


0.00 


7.8 


60.0 


0.130 


1.00 


1670 


M2 


King 


0.00 


8.7 


60.0 


0.145 


1.00 


1990 


KlOO-OOa 


Plum 


0.00 


6.6 


51.0 


0.129 


1.00 


1440 


KlOO-OOb 


Plum 


0.00 


6.6 


51.0 


0.129 


1.00 


1430 



Table 2. As for Table[T]but now showing properties of the models at T = 12 Gyr. Column 1 gives 
the label assigned to each model. Column 2 gives the half-mass radius followed by an estimate of 
the error in this value. This is repeated in columns 4 and 5 for the 2-dimensional projected half- 
light radius. The values in columns 2 and 4 are computed from the average of five snapshots taken 
over a 300 Myr period centred on 12 Gyr. Errors in columns 3 and 5 are based on the difference 
between the average and the greatest outlier. Columns 6-8 give the tidal radius, the mass that 
remains bound within the tidal radius and the half-mass relaxation timescale. All lengths are in 
pc units, masses in Mq and relaxation times are in Myr. 



Label 




error 


rh,l 


error 


rt 


M 




Nl 


10.5 


±0.1 


5.6 


±1.1 


106.0 


32 470 


3240 


N2 


17.3 


±0.1 


10.5 


±1.0 


103.0 


29 940 


6680 


N2b 


17.3 


±0.1 


11.1 


±1.8 


103.0 


29 960 


6400 


N3 


23.0 


±0.1 


15.4 


±1.0 


97.0 


25 410 


9440 


Ml 


8.1 


±0.2 


5.6 


±1.2 


38.4 


15 430 


1410 


M2 


5.6 


±0.1 


3.9 


±1.4 


32.5 


9 420 


620 


KlOO-OOa 


5.9 


±0.1 


1.0 


±1.7 


35.4 


15 780 


890 


KlOO-OOb 


6.5 


±0.1 


3.2 


±0.7 


35.2 


15 470 


1020 



(Spitzer 1987; Binney & Tremaine 1987). For a cluster to 
reach the dynamically-evolved state of core-collapse it has 
been shown that at least ten half-mass relaxation times must 
have elapsed (e.g. Baumgardt, Hut & Heggie 2002; Hurley 
2007) . The half-mass relaxation timescales corresponding to 
the initial states of the clusters are given in Table [1] for our 
models. In Table[5]the timescales for the clusters at an age of 
12 Gyr are given. We can see that of models Nl, N2 and N3 
it is only model Nl, the initially most compact cluster, that 
is likely to reach core-collapse within 20 Gyr of evolution. 
This is borne out in Figure [T] which shows the evolution of 
r c and rh for the three models. The tidal radius of model 
N3 is also shown for reference. Indeed, model Nl is close to 
core-collapse when the simulation is stopped at 20 Gyr. By 
comparison, the core radii of models N2 and N3 are showing 
no signs of decreasing to the small values associated with 
core-collapse (at 12 Gyr r c for N3 is the same size as rh for 
Nl). Based on the half-mass relaxation timescales such slow 
evolution for models N2 and N3 is expected as at 20 Gyr 
they have evolved for less than four and three half-mass 
relaxation times, respectively. 

It can be seen in Figure[T]that for all three models there 
is an initial rapid phase of increasing rh. This corresponds to 
the phase of violent relaxation associated with stellar evolu- 



tion mass-loss from massive stars in which there is an overall 
expansion of the cluster. Afterwards the evolution of rh pro- 
ceeds slowly and for models N2 and N3 rh is effectively con- 
stant over Gyr timescales. In model Nl we do see a continual 
increase in rh as the outer regions of the cluster expand in 
response to the shrinking core. At all times rh for N3 re- 
mains larger than for N2 which in turn remains larger than 
for Nl: the initial hierarchy of cluster size is maintained and 
the initially more extended cluster remains so. It is true that 
the spread in rh decreases with age - rh for N3 is initially a 
factor of three greater than for Nl compared to a factor of 
two greater at 20 Gyr - but even so the initial cluster size 
scale has an important bearing on the subsequent evolution. 

As described by Heggie & Hut (2003) an interesting 
way to understand the evolution of a star cluster is through 
ratios of the key radii (r c , rh and r t ). In Figure [2^ we re- 
produce Fig. 33.2 of Heggie & Hut (2003) and use it to give 
a basic overview of the expected behaviour (see Ch. 33 of 
Heggie & Hut 2003 for more background and detail). The 
evolution paths are distinct for clusters initially contained 
well within the tidal radius (path A in Fig. [2^) or filling the 
tidal radius (path B). For the former case all regions of the 
cluster expand during the initial period of violent relaxation 
so that r c /rh ~ constant while rh/r t increases (as r t changes 
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Figure 2. Evolution of star clusters in the r^/rt versus r c /r^ plane. Panel a) is for illustrative purposes (see text for an explanation of 
the labels) while panel b) shows Models 3 (open triangles) and 7 (grey crosses) of Hurley et al. (2004). In panel c) we show models Nl 
(crosses), N2 (open circles) and N3 (open triangles) of this work. Finally, panel d) shows models Ml (open circles) and M2 (solid stars) 
of this work as well as Models KlOO-OOa (open triangles) and KlOO-OOb (crosses) of Hurley (2007). 



only minimally). Thus the cluster evolves along path Al. For 
the tidally-filling cluster, where in the extreme we assume 
that rh follows the evolution of r t , we would instead ex- 
pect an increase in r c /rh while rh/r t remains constant (path 
Bl). Following violent relaxation the cluster then enters a 
long-lived period of dynamical evolution driven by two-body 
relaxation where mass-segregation develops and the cluster 
evolves towards core- collapse. During this time r c decreases 
markedly, rh changes only slightly by comparison and rt de- 
creases slowly. This gives paths A2 and B2 as both regimes 
head towards low r c /rh. Post-collapse evolution is character- 
ized by a central energy source (e.g. hard binaries) heating 



the core which will oscillate in size but with a time-average 
that varies little. For an isolated cluster we would expect 
rh to increase during this phase so any cluster that is not 
yet tidally-filling will move along path A3 to join B3. The 
tidally-filling clusters are expected to have little evolution 
in the r c /rh-rh/r t plane after core-collapse. We can readily 
verify this by looking at the evolution of some actual mod- 
els. First we start with two models of N = 30 000 stars from 
Hurley et al. (2004) where one model started well within r t 
and the other started with r max = r t . These are shown in 
Figure [2p and we clearly distinguish path A and path B-like 
evolution before both clusters end at a similar point. In real- 
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ity we would expect most clusters to evolve between the two 
extremes of paths A and B as they head towards B3, partic- 
ularly during violent relaxation where clusters will expand 
and can move along path CI to join B2. We see this to some 
extent when we look at the evolution of models Nl, N2 and 
N3 in Figure [2J;. Model N3 in particular first evolves across 
the phase space before starting down path B2. There is also 
a clear distinction between the evolution of the three mod- 
els, residing at well separated r^/rt values as they start to 
move down in r c /r%. Of course, owing to the large half-mass 
relaxation timescales, the models do not get the opportunity 
to evolve completely through the parameter space. 

We have now seen that it is possible to reach large 
values for clusters evolving in a weak NGC 6822-like tidal 
field and that quite distinct values can be obtained by 
clusters with different initial sizes. But can distinct rh val- 
ues be obtained by other means? The first possibility we 
explore is the inclusion of a primordial binary population 
(Model N2b). Figure [3] compares the rh evolution of models 
N2 and N2b which are identical in setup except for a 5 per 
cent primordial binary population in the latter. We see that 
the rh evolution is indistinguishable. This is also true for the 
evolution of bound cluster mass (see Tabled and other gen- 
eral quantities. We note that as an accuracy check models 
N2 and N2b were both performed twice with different initial 
random number seeds and the variation between different 
realisations of the same model was less than the difference 
between the two model types. For these extended clusters 
it is not surprising that the addition of binaries makes lit- 
tle difference to the evolution. The long relaxation times for 
Models N2 and N2b (see Tables [l] and [2]) mean that the col- 
lisions! extraction of binding energy from the binary orbits 
will not be efficient at heating the cluster. 

The second possibility for internal evolution creating 
differences between models relates to the formation of tight 
BH-BH binaries which then act as a central energy source 
to heat the cluster. Hurley (2007) showed that the forma- 
tion of one such long-lived BH-BH binary could double the 
r c /rh ratio compared to a similar model which did not form 
a BH-BH binary. However, we compare the rh evolution of 
these models (KlOO-OOa and KlOO-OOb from Hurley 2007) 
in Figure [3] (and in Tables [5} and we see no clear distinc- 
tion. Mackey et al. (2008) took this further and contrasted 
the evolution of clusters with no BH-BH binaries to that of 
clusters which retained a large number of post-supernovae 
BHs (~ 200) that subsequently sank to the cluster centre 
and formed BH-BH binaries (as many as five such binaries 
present at any one time). The focus was on star clusters in 
the Large Magellanic Cloud and as such a tidal field the same 
as for our NGC 6822 case was used (but with i? gc = 6 kpc 
rather than 10 kpc). They found that the inclusion of the 
BH-BH binaries could increase rh by as much as a factor of 
two by the time that the model with no BH-BH binaries had 
reached core-collapse (compared to a corresponding factor 
of 20 increase in r c ). The rh behaviour first started to di- 
verge after 1 — 2 Gyr of evolution, corresponding to roughly 
one half-mass relaxation time, with the expansion driven 
on the shorter relaxation timescale of the centralised BH 
population. However, we note that the most extended of 
these models were in the advanced stages of dissolution at 
a Hubble time. Another possibility is one that has gath- 
ered much attention of late, namely the question of whether 
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Figure 3. Comparison of the half-mass radius evolution for mod- 
els N2 and N2b (solid lines), Ml and M2 (dashed lines), and 
KlOO-OOa and KlOO-OOb (dotted lines) of Hurley (2007). 



or not some star clusters harbour intermediate-mass black 
holes (IMBHs). Gill et al. (2008) compare the rh evolution 
of models with and without an IMBH and find no significant 
difference until well after core-collapse and even at very late 
times the difference is still less than a factor of two. Baum- 
gardt, Makino & Hut (2005) looked at the effect of increas- 
ing IMBH mass and found an increase in rh of 15 per cent 
at most. It should be noted that the maximum black hole 
mass included in the models to date is 1 000 Mq and that 
the heating produced by significantly more massive IMBHs, 
if indeed they exist, is yet to be documented. 

We next look at the evolution of our models Ml and 
M2 which were evolved in the stronger M31-like tidal field. 
These models started with different initial density profiles so 
provide an opportunity to look at how the choice of a Plum- 
mer or King profile affects the rh evolution. This is shown 
in Figure [3] (also in Figure [4} and we see that at various 
stages in the evolution the difference can be up to 50 per 
cent. The radius evolution of models Ml and M2 is studied 
in more detail in Figure [4] Comparing this to Figure [1] we 
clearly see that the stronger tidal field drives more rapid 
evolution for the M31 models relative to their NGC 6822 
counterparts. Indeed, both Ml and M2 reach core-collapse 
prior to 20 Gyr. The model with the King density profile 
evolves more rapidly. This is primarily owing to a greater 
central density of stars in the initial model which led to a 
greater rate of dynamical interactions, more mass lost across 
the tidal boundary in the early stages, and consequently a 
reduced relaxation timescale. However, in the M-rh plane 
(Figure [SJ) the evolution of models Ml and M2 is very sim- 
ilar. Interestingly both Ml and M2 appear to keep rh ap- 
proximately constant after core-collapse, behaviour that was 
noted by Kiipper, Kroupa & Baumgardt (2008) in their 
models. We also see that models Ml and M2, as well as 
models KlOO-OOa and KlOO-OOb of Hurley (2007) all evolve 
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Figure 4. Evolution of the half- mass radius, rh, and the core- 
radius, r c , for the M31 models Ml (solid lines) and M2 (dashed 
lines). The tidal radius of model Ml is also shown. 
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Figure 5. Cluster mass versus half-mass radius for the models 
in this work as well as models KlOO-OOa and KlOO-OOb of Hurley 
(2007). For each evolution track the point at 12Gyr is denoted 
by an open circle. 



towards a similar value of rh — 5 pc even though they start 
with distinct initial conditions or follow different evolution 
paths in other respects. The evolution in the r c /rh-rh/r t 
plane is also shown for these models in Figure [2}! where 
they all lead to a final r^/r t — 0.18 and appear to evolve 
intermediate to the N2 and N3 cases (although for N3 it is 
not dynamically evolved enough to determine if the path is 
distinct from Ml or not). 

Figure [S] looks at the evolution tracks of our model clus- 
ters in the M-rh plane. We see that for any particular mass 
the initially extended clusters in the weaker tidal field have 
the greater half-mass radii. Therefore, clusters at the same 
R gc in M31 and NGC6822 will follow distinct paths in the 
M^rh plane, assuming they were all close to filling their 
tidal radii initially. This is to be expected. However, the sit- 
uation changes if we look at the cluster mass versus rh/r t 
ratio as in Figure [6] Now the tidally-filling clusters in M31 
and NGC 6822 follow very similar evolution paths. Kiipper, 
Kroupa & Baumgardt (2008) showed that clusters evolve to- 
wards a common post-collapse sequence in this plane. Our 
models would appear to back this finding (a point we will 
return to in Section [4j and the sequence in our case corre- 
sponds to a constant or slightly increasing rh/r t ratio. For 
models Nl and N2 it is too early to tell if they too will join 
the sequence when they eventually reach the post-collapse 
regime but indications are that this is likely. 

In Figure [7] we look at how the bound cluster mass 
changes with cluster age for the various models in this work. 
We also include the KlOO-OOa and KlOO-OOb models of Hur- 
ley (2007) for comparison. As we would expect the clusters 
in the weaker NGC6822-like tidal field lose stars/mass at 
a slower rate compared to clusters orbiting within a galaxy 
that is ten times more massive. Also, for models Ml and M2 
(where M2 is the lower of the dashed lines) we once again 
see that evolution is more rapid for the Wo = 7 King model 
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Figure 6. Cluster mass versus the ratio of the half-mass to tidal 
radius for the models in this work as well as models KlOO-OOa 
and KlOO-OOb of Hurley (2007). 

compared to the Plummer model, although the slopes are 
similar from about 5 Gyr onwards. For the NGC 6822 model 
clusters we see that more mass is lost by any particular 
time if the cluster was initially closer to filling its tidal ra- 
dius (noting that Nl is the upper of the solid lines). This 
corresponds to the finding of Gieles & Baumgardt (2008) 
that clusters that are initially well within the tidal radius 
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Figure 7. Cluster mass (scaled by the initial mass) as a function 
of age for the models in this work as well as models KlOO-OOa and 
KlOO-OOb of Hurley (2007). 




Figure 8. The escape rate (k e : see text for details) as a function 
of cluster age (scaled by the half-mass relaxation time) for the 
models in this work as well as models KlOO-OOa and KlOO-OOb of 
Hurley (2007). 



take a much greater number of half-mass relaxation times 
to lose half of their mass than clusters that fill their tidal 
radius. Following Giersz & Heggie (1997) an expression for 
the mass-loss rate can be written as: 



where In A is the Coulomb logarithm and we take A = 0AN. 
For their A-body models of N = 500 stars Giersz & Heggie 
(1997) found a value of 1.3 for fc c . In Figure [8]we plot k c as 
a function of the number of half-mass relaxation times that 
have elapsed. We see that the behaviour for all of the clusters 
that initially fill their tidal radii is similar, in particular N3 
follows Ml, and that fc e decreases as the clusters become 
more evolved. We also see that fc e is lower for clusters with 
smaller initial rh/r t , a result that is predicted analytically 
by Gieles & Baumgardt (2008) and previously demonstrated 
by the Monte Carlo models of Spitzer & Chevalier (1973). 
In this work we use Eq. [3] and Figure [8] simply as a way 
of comparing the mass-loss evolution of the various models. 
A detailed examination of how dissolution times for star 
clusters in tidal fields can be expected to scale with cluster 
parameters can be found in Baumgardt & Makino (2003) 
and Gieles & Baumgardt (2008). 

We also remind the reader at this stage that the values 
of rh are based on three-dimensional data. If instead we cal- 
culated the half-mass radii from a two-dimensional projec- 
tion a reduction of about 25 per cent would result (Fleck et 
al. 2006). Furthermore, Hurley (2007) showed that projected 
half-light radii could be as much as half that of the corre- 
sponding half-mass radii for dynamically evolved clusters. 
In Table [2] we show both the half-light (two dimensional) 
and half-mass (three-dimensional) radii for the models at 
an age of 12 Gyr. The errors in both values are also shown. 
Most obvious is that the errors in rh,i are greater than in 



7"h reflecting that the former quantity fluctuates on a short 
timescale resulting from a sensitivity to the movement of 
massive stars in combination with events such as the evolu- 
tion of bright giants to become dim white dwarfs. We can 
find particular times when 7~h,i is as low as 50% that of rh. 
However, on average the reduction is in the 30-35% range. 
It might also be expected that the r^j/rh ratio would be 
smaller for more evolved clusters with a greater degree of 
mass-segregation (say Ml compared to N3) but overall we 
see only a weak trend for this. 



4 DISCUSSION 

We have shown with our model N3 that clusters in the outer 
regions of a galaxy such as NGC 6822 can naturally evolve 
to have large half-mass radii provided they begin by filling 
their tidal radii. The half-light radii in projection can ex- 
ceed 15 pc which is within the range observed for extended 
clusters. Furthermore, extended clusters are observed out to 
~ Ukpc in projection from the centre of NGC 6822 which 
corresponds to a three-dimensional i? gc of ~ 15 kpc. Our 
models were evolved at R gc = 10 kpc so following equation[T] 
we can scale the tidal radius by a factor of 1.5 and poten- 
tially the half-mass and half-light radii by a similar factor. 
Thus half-light radii in excess of 20 pc are possible. Also con- 
sider that our models started with N = 100 000 and have a 
mass of about 30 000 Mq remaining at 12 Gyr. If some of the 
actual clusters are in fact more massive, say twice as mas- 
sive, then following equation [T] again gives another potential 
25% increase in the quoted half-light radii. Therefore, for 
a dwarf irregular galaxy such as NGC 6822 extended clus- 
ters with projected half-light radii in the range 10 — 30 pc 
can be expected to naturally occur (for a range of masses 
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and orbits). Such clusters will evolve slowly. They will be 
dynamically young, having evolved for approximately two 
half-mass relaxation times, and will not show any signs of 
core- collapse. 

Any clusters in the outer regions of NGC 6822 with 
half-light radii less than 10 pc will need to have started 
life relatively compact. Certainly they would not be ex- 
pected to have filled their tidal radii initially. For a cluster 
in NGC 6822 to exhibit core-collapse within a Hubble time 
an initial filling factor of less than 30% is required, or an 
orbit at R sc much less than 10 kpc. We have shown that for 
a fixed i? gc differences in initial sizes can translate to dif- 
ferences in half-light radii at 12 Gyr by the same factor. We 
have also shown that primordial binaries do not affect the 
evolution of the half-mass and half-light radii. This is not 
surprising for dynamically young clusters and has also been 
demonstrated for more evolved clusters in stronger tidal 
fields (Hurley 2007). The internal evolution characteristic 
that has the most potential to impact the half-light radius 
evolution is the formation of a central population of BH-BH 
binaries. Expansion of the half-light radius by this mecha- 
nism was shown by Mackey et al. (2008) to start taking effect 
after one half-mass relaxation time, with the effect increas- 
ing to a maximum expansion factor of two by core-collapse 
(comparing clusters with and without BH-BH binaries). So 
this could be a mild factor for the most extended clusters 
in NGC 6822 and potentially play a larger role in determin- 
ing the observed half-light radii of initially more compact 
clusters. However, this relies on a source of BHs being re- 
tained in the clusters. Overall we find that differences in 
initial cluster size are the most likely explanation for any 
measured differences in half-light radii for clusters that oth- 
erwise appear similar. The stochastic effect of the presence, 
or not, of a central population of BH-BH binaries can then 
possibly play an added role. 

Our question of "how do similar clusters obtain differ- 
ent structure?" now seemingly becomes "can clusters be ex- 
pected to be born with different sizes?". Baumgardt et al. 
(2010) asked a similar question after showing that Milky 
Way GCs beyond R gc of 8 kpc fall into distinct groups of 
compact, tidally-underfilling clusters and more extended, 
tidally-filling clusters. They concluded it was more likely 
that the second group formed with large half-mass radii 
rather than evolving from a more compact state. Examining 
the size distribution of globular clusters in a variety of Local 
Group dwarf galaxies, Da Costa et al. (2009) suggest that 
there may be two distinct modes of cluster formation with 
typical half-light radii of 3 and 10 pc, analogous to the find- 
ings of Baumgardt et al. (2010). It has been suggested by 
Ballesteros-Paredes et al. (2009) that tidal forces may play 
a role in determining the rate at which stars form within 
molecular clouds, through compression and possible disrup- 
tion of the clouds. Their calculations show that the resultant 
star formation efficiency of a particular cloud depends on its 
position and orientation in the host galaxy. This model is 
put forward as an alternative to the theory that magnetic 
field strength and radiative feedback regulate the degree of 
star formation within a cloud (Price & Bate 2009). Either 
way, if differences in star formation efficiency can exist be- 
tween proto-clusters forming in the cores of molecular clouds 
then Baumgardt & Kroupa (2007) have shown that differ- 
ent 9"h/r t values can result for the bound clusters remaining 



after residual gas is expelled. Furthermore, simulations con- 
ducted by Elmegreen (2008) have demonstrated that diffuse 
star clusters can form directly in regions where the back- 
ground tidal forces are low. Here star formation efficiency 
is discussed in terms of turbulence (or Mach number) and 
cloud pressure (or column density). Different combinations 
produce clusters of different densities with high-pressure re- 
gions making the densest clusters while clouds with low den- 
sity and high Mach numbers give rise to diffuse clusters. 
The latter are more likely to form in the outer regions of a 
galaxy where low-density clouds are expected to reside and 
the diffuse clusters will not be quickly disrupted (Elmegreen 
2008). So there are possibilities for variations in environ- 
ment and /or the properties of molecular clouds producing 
different star formation efficiencies and a range of r^/rt val- 
ues for young star clusters. What is not immediately clear 
is whether or not clusters forming in the same environment 
can have a range of r^/rt values owing to variations in cloud 
properties (size, shape, orientation) alone. 

We now turn our focus to the M31 extended cluster 
system. Our model Ml which started tidally-filling demon- 
strates the maximum size that a cluster of M ~ 15 000 Mq 
orbiting at R sc — 10 kpc can be expected to have at an 
age of 12 Gyr. This is not within the extended cluster size 
range. Of course we are not claiming that our model Ml 
is a direct representation of a typical star cluster in M31. 
For a start, the mass of the model is likely to be an un- 
derestimate compared to actual M31 clusters. The extended 
clusters observed in M31 have masses that range from ~ 10 4 
to ~ 10 Mq (Huxor et al. 2010a: note that masses for ECs 
in NGC 6822 are not yet known but are expected to cover a 
similar range). Thus we have direct models at the lower end 
of this range but could be as much as a factor of ten too small 
at the upper end. A factor of ten increase in mass will lead 
to a larger tidal radius, but only by a factor of two and if the 
internal radii scale in the same way this gives a projected 
half-light radius in the range of 10 pc (see Table [2}. This 
is right at the low-end of the extended cluster size range. 
However, this still does not help when we are looking at two 
clusters of the same mass (or luminosity) at the same R gc 
and with very different half-mass radii, as is observed in the 
inner regions of M31 (Huxor et al. 2010a). Thus it is hard to 
believe that extended clusters observed in the inner regions 
of M31, which have half-light radii of about 30 pc, evolved 
in-situ. It is possible that extended clusters observed in the 
inner regions are actually further out and only appear in 
the centre owing to projection effects. Such an effect is not 
noticeable from the horizontal-branch levels of the few ex- 
tended clusters that have good colour-magnitude diagrams 
(Mackey et al. 2006) and if there is a projection effect push- 
ing these clusters into the foreground there is no reason why 
it shouldn't similarly apply to the compact clusters as well. 
On the other hand we expect that compact clusters will be 
more likely to survive in the inner parts of galaxies than ex- 
tended clusters and radial distributions of GCs in galaxies 
show that numbers increase towards the centre (see Brodie 
& Strader 2006, for example). So for compact clusters at 
least it is reasonable to believe that there is a sizeable popu- 
lation residing in the inner regions of galaxies such as M31, 
regardless of projection effects. 

It is also true that for our models we have used an 
i? gc that for M31 represents only the very inner edge of the 
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extended cluster locations. As discussed above, and shown 
by equation U in a smooth tidal field the tidal radius of 
a cluster scales linearly with _R gc . So the maximum possi- 
ble rh will also increase for clusters observed further out 
than lOkpc. The M31 extended clusters are observed out 
to R sc of 100 kpc and beyond which means a potential fac- 
tor of ten increase in size for clusters in the outer regions. 
We note though that other factors such as the relaxation 
timescale will be affected by a change in cluster size so it 
is not necessarily a simple case of scaling up a model value 
at a particular age. To gain some insight we can turn to 
Mackey & van den Bergh (2005) who have indeed shown a 
clear trend of increasing with increasing R gc for Milky 
Way globular clusters (see their figure 9). However, the re- 
lation is not linear: going from R gc of 10 to 100 kpc gives a 
factor of 2 — 3 increase in the typical cluster size (as mea- 
sured by rh). If we assume a similar scaling for M31 then 
we can expect to find clusters with rh,i ~ 30 pc in the outer 
regions. It is particularly interesting to note that the rh,i 
values derived for the M31 extended clusters (Huxor et al. 
2010a) show no obvious correlation with R gc . In combina- 
tion with our model results this infers that the inner M31 
extended clusters either resided at more distant locations 
for a large part of their lives before moving to their current 
locations or they were likely accreted as part of one or more 
dwarf galaxies (perhaps similar in nature to NGC 6822 or 
the dwarf elliptical Scl-dEl). This latter scenario fits with 
the picture sketched by Da Costa et al. (2009) where the 
formation of extended clusters can be expected to be more 
common in dwarf galaxies where tidal effects are less dis- 
ruptive and that extended clusters seen in larger galaxies, 
M31 for example, were accreted in dwarfs that subsequently 
disrupted. Adding to this argument is evidence of a genuine 
physical association between the extended clusters (as well 
as most compact clusters) and coherent stellar streams in 
the halo of M31 (Mackey et al. 2010). 

Our models so far have been evolved within a smooth 
background potential represented by a point-mass galaxy. 
This simplification may cause discrepancies when comparing 
our models to clusters that evolve in host galaxies with dis- 
tinct structures. In particular, both the Milky Way and M31 
exhibit pronounced bulge, disk and halo components which 
means that the effects of bulge and disk shocking need to 
be considered. Gnedin, Lee & Ostriker (1999) have demon- 
strated with Fokker-Planck simulations that gravitational 
shocks caused by the time-varying tidal forces as clusters 
pass near the central Galactic bulge or through the Galac- 
tic disk can significantly reduce core-collapse and dissolution 
times. The cluster binding energy is reduced by the shocks so 
that diffuse clusters are the most susceptible to disruption. 
However, the effect of the shocks is strongly dependent on 
the Galactocentric radius. Bulge shocking dominates clus- 
ter evolution for R gc of 2 kpc or less (Gnedin & Ostriker 
1997) but drops off quickly away from the central nucleus, 
while disk shocking has a minimal effect for _R gc beyond 
about 8 kpc (Vesperini & Heggie 1997) owing to lower disk 
density and less frequent passages. To check this we have 
repeated model Ml using a three-dimensional Galactic po- 
tential: a point-mass bulge, a Miyamoto & Nagai (1975) disk 
(with a = 4 kpc and b — 0.5 kpc) and a logarithmic halo po- 
tential (see Aarseth 2003 and Praagman, Hurley & Power 
2010 for the NB0DY6 implementation). Masses of 5 x 10 10 and 



1.5 x 10 Mq were chosen for the bulge and disk, respec- 
tively (Xue et al. 2008), with the halo mass chosen so that 
the mass interior to 10 kpc was the same as for model Ml. 
The cluster was evolved at 7? gc = 10 kpc perpendicular to 
the disk and there were no obvious changes in evolution of 
cluster mass or the structural parameters compared to Ml. 
Thus we are confident that our use of a point-mass galaxy 
has little bearing on our results for clusters evolving at R gc 
of 10 kpc or greater. Conversely, as suggested by Da Costa 
et al. (2009), the action of tidal shocks makes it unlikely 
that extended clusters would be found at smaller galactic 
radii in large disk galaxies such as M31. A full investigation 
of the effect of three-dimensional galaxy potentials on the 
evolution of star clusters will be the focus of future work, 
noting that the scale-length of the disk of M31 (~ 6 kpc) 
is larger than for the Milky Way disk (~ 4 kpc: Yin et al. 
2009). 

Another avenue for future study is the effect of the 
initial cluster density profile on the subsequent evolution. 
Trenti, Vesperini & Pasquato (2010) find that King models 
with different initial concentrations evolve to a similar struc- 
tural state within a strong tidal field. However, the differ- 
ences we observed between our models Ml and M2 (evolved 
in the M31-like tidal field) means that the effect of chang- 
ing the initial core-density should be investigated for models 
evolving in a weaker NGC 6822-like tidal field. 

We should also note that the results and conclusions we 
draw regarding extended clusters are based on the assump- 
tion that these are essentially normal star clusters where the 
dynamical processes are not influenced by the presence of 
any significant dark matter component. If extended clusters 
are found to contain dark matter then our results to date are 
not applicable. The only direct constraints on this possibility 
thus far are the measurements for one M31 extended cluster 
reported by Collins et al. (2009). Their derived mass-to- light 
ratio of M/L — 6.7lg 5 7 Mq/Lq rules out very large ratios 
as seen in most dwarf spheroidals. As stated by Collins et 
al. (2009) the value is most consistent with a typical star 
cluster. However, the possibility of M/L as high as 20 still 
exists so the issue of whether or not extended clusters con- 
tain dark matter is yet to be resolved. The case against dark 
matter is strengthened by Jordi et al. (2009) who measured 
radial velocities of red giants in Palomar 14, a diffuse GC in 
the outer regions of the Milky Way, to derive a modest ratio 
of M/L ~ 2. This is compatible with the value expected for 
a normal stellar initial mass function with no dark matter. 

Finally, our models build on a number of important re- 
sults from previous work regarding general star cluster evo- 
lution. A prime example is the work of Kiipper, Kroupa 
& Baumgardt (2008) who used TV-body models to demon- 
strate that post-collapse clusters asymptotically evolve on 
to a common sequence where the Vh/rt ratio depends only 
on the mass remaining in the cluster and not on the initial 
conditions. They considered open clusters starting with up 
to N ~ 32 000 stars and different binary fractions, density 
distributions, and R sc values. Our models extend this 
work to include a distribution of stellar masses and stel- 
lar evolution for N — 100 000. In doing so we confirm the 
findings of Kiipper, Kroupa & Baumgardt (2008). The com- 
mon sequence shows a rh/r t ratio that starts in the range 
of 0.15 — 0.2 and increases slightly with decreasing mass. 
Related to this the same authors found that clusters in a 
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tidal field show an equilibrium half-mass radius after core- 
collapsc. This was 2 pc for their standard set of models, not- 
ing that the particular value will increase for higher initial 
mass and larger i? gc . Our models Ml and M2 have a similar 
tidal tidal field but larger initial mass. These models also 
show an equilibrium half-mass radius and as expected it is 
larger, ~ 5 pc, with the inclusion of stellar evolution likely 
contributing to some of the increase. This corresponds to 
a projected half-light radius of about 3pc which fits with 
the typical value observed for compact GCs in the Milky 
Way (Baumgardt et al. 2010) and M31, and indeed almost 
all galaxies where we see GCs (Da Costa et al. 2009, for 
example). Of our models in the weaker tidal field only Nl 
gets close to core-collapse and indeed has a larger half-mass 
radius at this point, ~ 10 pc, but in this case appears to still 
be rising. 



5 CONCLUSIONS 

Using a direct JV-body code we have followed the evolution 
of star clusters with different initial sizes in a tidal field ap- 
propriate for a galaxy such as the dwarf irregular NGC 6822 
or the Large Magellanic Cloud. We also looked at the effect 
of increasing the galaxy mass by a factor of ten, appropriate 
for larger galaxies such as M31 and the Milky Way, on the 
evolution of clusters which initially fill their tidal radii. Our 
main findings can be summarised as: 

• extended clusters with projected half-light radii of up to 
30 pc in galaxies such as NGC 6822 can result from standard 
star cluster evolution provided that the clusters are close to 
filling their tidal radii initially; 

• internal evolution factors such as the proportion of pri- 
mordial binaries or the formation of BH-BH binaries do not 
produce enough expansion so formation as an extended ob- 
ject followed by standard cluster evolution is the most likely 
path for the extended clusters recently observed; 

• observations of clusters of similar mass and age but sig- 
nificantly different half-light radii most likely indicate that 
the clusters started life filling their tidal radii by different 
factors; 

• in a galaxy with a strong tidal field, such as M31, the 
maximum half-light radius expected for clusters that evolve 
at a galactocentric radius of lOkpc is of the order of 10 pc; 

• extended clusters could form naturally in the outer re- 
gions of M31 (at ~ 100 kpc) but any extended clusters ob- 
served in the inner regions likely formed in one or more dwarf 
galaxies that were subsequently accreted by M31. 
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